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ABSTRACT 


This is a final study of a radiational model for use in the FNWC 
prediction system. This model utilizes original FNWC data over oceanic 
gridpoints which were adapted to the FNWC five-layer initial data 
analysis in four mid-seasonal studies for the data year 16 January 
1974, 16 April 1974, 16 July 1974 and 16 October 1973. 

The radiational model parameterizes the large-scale cloud amounts 
in two layers. The primary objective was to tune model albedo values 
to those taken from satellite climatology by changing the cloud-reflec- 
tance coefficients used in previous studies. The albedo values were 
tuned for least squares deviation relative to satellite climatology 
albedos. The present model using tuned albedos better verified the 
radiational balance at the top of the atmosphere against satellite 
climatology than did pre-existing untuned models, both on a seasonal 


aaa annual basis. 
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Zs M@NTRODUCTION 


This thesis is a refinement of previous radiational models described 
Pein) spacth (1975), (B) Meyers (1975), (C) Beahan (1975) and (D) 
Warner (1974) for use in the Fleet Numerical Weather Central (FNWC) pre- 
diction system. This study has as a primary objective, the comprehensive 
re-examination of the radiational physics in layers comprising the ocean- 
atmosphere system. FNWC atmospheric soundings defined at constant pres- 
Sure levels for the four mid-seasonal dates of the "data year" 1973-74 
aS previously examined in studies (A,B,C,D) were utilized in this study. 

The application of the FNWC radiative model may be made to any scale 
of analysis for which there is adequate resolution of the temperature 
and moisture data in the vertical. In the horizontal, the reliability 
of the data used here is consistent with that of the FNWC interpolation 
to gridpoints in the analysis procedure. Temperature and dew-point data 
in radiative soundings are typically reported to the nearest tenth of a 
degree. The radiative computations made here are applied to FNWC grid- 
points and are designed to make a one-hour forward-time step at gridpoints 
in the FNWC primitive equation forecast model, with special adaptions to 
O-levels. 

The specification of cloud amounts in two designated layers, one at 
a mid-level and the other at a low-level, has an important influence on 
the radiative-model dispositions, both in the short- and long-wave 
Spectral regions. Initially in (A,B,C,D), the specification of the 
fractional amounts of CL. and CL. had been based on large-scale formula- 


1 2 


tions developed by Smagorinsky (1960), but during the course of these 


ie 





studies it was found more realistic to modify initial CL-values to 

Ci = 273 CL. This reduction in cloud amount was also used here as it 
had been chosen to prevent CL' from exceeding unity and to afford 
better agreement with satellite climatology, such as planetary albedo, 
for compatible data periods. 

With the reduced cloud coverages CL, ' and CL.' from the earlier 
studies (A,B,C,D), it was possible to obtain reasonably close agreement 
in the computed terrestrial net flux at the top of the atmosphere and 
that observed for comparable NIMBUS III subsatellite points and data 
periods (Raschke, Von der Haar, Bandeen and Pasternak, 1973). However 
even with the reduced cloud coverages CL! and cL." the computed 
Planetary albedo remained generally excessive, particularly in tropical 
latitudes. Hence it became a major objective of this particular study 
to modify empirically the reflective capability of the cloud layers. 
This was done by defining a general factor f so that the initial choices 


(after Rodgers, 1967) of cloud reflectances R were modified to R‘' where 
R' = fR, 


Systematic substitution of the cloud reflectances R, wherever they 
entered the solar disposition equations, by R' then led to a relation- 
ship between the global albedo and f. Utilization of the least squares 
technique to minimize the differences between satellite and model albedos 
Over a geographic sample of points led to best-fit value of the "tuning- 
factor" £. Separate values of f£ were deduced by least squares for each 
season, and subselections were deduced for the tropical and extratropical 


areas, respectively. 
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The modified solar cloud-reflectances improved the agreement 
between satellite and radiative model albedos in both geographic 
areas insofar as net incoming insolation was concerned. The terres- 
trial net flux at the top of the atmosphere was unaffected by the 
choice of f, while the use of CL! and CL! as specified gave good 


agreement with satellite terrestrial net flux data over the geo- 


graphic range and the time-scales concerned. 
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Ii. DATA PREPARATION 


eens FIELDS 
1. General Considerations 

The initial temperature and humidity data used in this study 
were arranged in the form of soundings taken along four oceanic meri- 
dians (Fig. 1) of the Fleet Numerical Weather Central (FNWC) Northern 
Hemisphere mid~-seasonal analyses for 16 October 1973, 16 January 1974, 
16 April 1974 and 16 July 1974. Oceanic locations for these computa~ 
tions were chosen because: 

(a) constant O~-surfaces (where O = =) of the FNWC primitive 
equation system are nearly identical constant pressure levels. 

(b) The maritime-area soundings are more likely to be system- 
atically representative of the set of zonally-distributed gridpoints 
than over land. 

The three meridians selected over the Pacific Ocean were lo- 
cated at 125W (25 soundings), 170W (25 soundings) and 145E (17 sound- 
ings). The Atlantic Ocean meridian was located at 35W (26 soundings). 
This method of selecting "soundings" along the indicated meridians 
of the FNWC polar stereographic map made it unnecessary to employ spa- 
tial interpolation between original data gridpoints along the meri- 
dians. Data along line 3 in the Pacific was not extended southward of 
gridpoint (9,55) because they fell over land masses (New Guinea and 
Northern Australia) where the surface temperatures and other sounding 


features were unrepresentative of the oceanic values. 
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‘Figure 1. FNWC polar stereographic grid and meridians (lines l, 2, 


3, and 4) selected for study. The longitudes A are shown 


for each meridian as well as the extent considered of each 
"meridian. 

















ee Daca LTheeacment 
a. Original Soundings and Modifications 

The gridpoint soundings were taken from the original FNWC 
63-by-63 Northern Hemisphere analyses of T(p) and of T - Ty (the dew 
point depression) at standard pressure levels up to and including 
p ~100 mb. Examples of such original soundings were shown as Table 
I(a) in Meyers (1975, p. 24). Subsequently each original FNWC sound- 
ing was transformed --in previous studies (A, B, C, D) of this series 
--into what has been termed the radiative sounding having the format 
shown in Table I. The data levels of the radiative sounding contains 
essentially the five FNWC predictive O-levels (dotted levels in Fig. 2). 

At each gridpoint selected, the original FNWC humidity 
soundings were given in the form of five dew point depressions over 
the analyses levels from 925 mb to 400 mb. At the surface (level 
k = 10) the standard instrument level vapor pressure, Cair’ was trans- 


formed into the surface mixing ratio, by means of 


Rays 


aay) C2lcu ates, / 100m. (a0 


To obtain radiative soundings as in Table I, it is neces- 


Sary to have water vapor and CO. absorber masses at certain required 


2 


k-level boundaries (Fig. 2). All radiative soundings in this study 


start at sea level with the approximation of surface pressure 


T= 1000 mb. Therefore, the eleven k-levels correspond closely to the 


FNWC levels P. =O oO Ore mow slo. . 592200. , LOO., O.0 mb and in 


Puden: co O) — Oren. is cy, OL, O20. 
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i elLeele Example Of a radiative sounding at gridpoint (1,1) for 
16 April 1974 with mixing ratio listed at odd k-levels (Fig. 2). 
Additionally, water-vapor and CO. absorber masses are also listed 
as these parameters have been modeled in the radiative theory pre- 
sented in this study. 


Mixing Absorber Masses 
Pressure Temp Ratio Water Vapor CO> 5 
(mb) Ge (g/kg) (gm/cm) (cm/cm*) 
1000 25.60 Uy aD 
900 16.09 NA 073} 
800 Tle 25.20 AD 53 
700 56 18)0) ois al 
600 2510 S623 BS 555 
500 aI 8) 2 
400 235.20 eS S Birr ss, 13.35 
300 253) (010 0.65 
200 =asele, Salen IRS Rei S12 
100 -80.70 0.03 SSP op 133.07 
O = 010 70, BO. a3 e355 


ibs, 





Level - 





identifers U(0), C60) O-level 
0 o=0 
Layer (0.2) Ieee eee [OP ae SS Se 0.1 
Qa 6=0.2 
(2,4) ei ae Ma a 0.3 


en Ke) Gh 


(4.6) 





(6,8) [poe Se aoe Se Iiqr-- 07 
0508 

(8,10) th hl aia a aaa 

i ig 10 


Figure 2, Five-layer radiative sounding used in this study. Levels 

are identified by their values on the k-scale, while layers are identi- 
fied by their level boundary indices in parentheses, e.g. (8,10). 
Pressure~-scaled water vapor and CO. mass increments AU and AC, respec- 
tively are integrated relative to the surface and the resulting U and C 
are carried at even levels. The temperature T is retained at all levels. 
Amounts of clouds CL; and CL 9 #£=in the layers shown have been para- 
meterized for consideration of their radiative effects. 
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b. Radiative Temperature Profiles 
mien Grrdeoint temperatures were listed at each mandatory 
mvweiwor Table I between 1000.,..., 100 mb (1.e., between k = 10,..., 
1). The temperature was assumed to be isothermal from 100 mb to 
0.0 mb. The temperature at level k = 10 was set equal to the FNWC 
listed sea-surface temperature. The radiative sounding temperatures 
for the remaining k-levels were obtained from either their correspond- 
ing listed temperature-level or by a three-point Lagrangian interpola- 
tion scheme [Egq. (2-1), Spaeth, 1975], to level k when the listed 
FNWC temperature profile did not include the value Ty 
c. Radiative Moisture Profiles 
Similarly, the moisture profiles of Table I have been ob- 
tained by an interpolative procedure over the original-level FNWC mix- 
ing ratios to those required at k-levels in a manner analogous to that 
discussed by Spaeth (1975, pp. 29-31). 
d. Pressure-Scaled Absorber Masses 
The pressure-scaled water vapor and the carbon-dioxide 
scaled absorber masses were calculated for the six even numbered k- 
levels (Fig. 2) using Eqs. (2-8, 2-9, 2-10 and 2-11) respectively as 
outlined by Spaeth (1975, pp. 31-33). These equations use essentially 


the mixing ratios of water vapor and of CO. at odd k-levels. 


B. CLOUD PARAMETERIZATION 

The relative humidities (RH) at levels k = 5 and k = 9 are used 
in the calculations of the fractional cloud covers CL, and CL. a 
layers (4,6) and (8,9) respectively (Fig. 2). The equations for 


Px I 





parameterization of the two fractional cloud amounts are as follows: 


(Gir eyo Ne.0 RH (S)) - 0.7) (2—2a) 


CL Zyemo.oS(kH(9)) = 2.0] (2-25) 


The bracketed part of the equations (after Smagorinsky, 1960) were 
Peauced by the 1/3 factor in an attempt to tune cloud amounts to obtain 
albedo values in closer agreement with the recent satellite radiational 
climatology of Raschke et al., (1973). The "2/3 CL" parameterization 
set forth in Eq. (2-2) is used here for estimating large scale radia- 
PEenelescttects only. hus small-scale convective activity and a priori 
climatological effects were not considered in specifying the form of 
Eq. (2-2). Tuning of the model-albedo values of this study by varying 
cloud reflectance coefficients with respect to Season and latitude will 


bemadiscussed in Section IV. 


C. CLOUD-AREA COVERAGES 


Fractional cloud amounts, CL. and CL, were computed at each grid- 


1 
point for levels k = 5 and k = 9 by Eqs. (2-2a,b). In addition the 


gridpoint area may be thought of as broken into random fractional seg- 


Mencs Of size 


W(0,0) = (1-CL, ) Cli Gina) (2-3a) 


z 


wherein there is a combination of clear-over-clear segments in the 
layers. Similarly, the gridpoint area has the fractional area of cloud 


coverage 


W(1,1) = CL, *CL. (oe) 


Ze 








of an upper overcast amount overlying a lower overcast. Likewise the 
area-combinations of overcast over clear, and clear over overcast 


areas in two layers, may be visualized as occurring with the weights 
W(1,0) = cL, *(1-CL,) (2-3c) 
and 


W(0,1) (2-3d) 


— * 
(1 CL) cL, 


respectively. 

For radiational computations it was useful to carry the relative 
weights or fractions of the gridpoint area exposed to the specified 
Cloud-layer combinations. Henceforth, the symbols denoted by W(0,0), 
fee), W(l,O) or W(0,1) indicate the fractionally overcast (1) or 
clear (0) cloud-area combinations in the indicated layers (Fig. 2), 


with the first index 1 or O referring to layer CL k = 5, and the 


1 
second to CL, k = 9, 

The usefulness of the cloud~area weighting device will be clarified 
in Sections III and IV, where the procedures for the terrestrial and 
solar radiational computations are discussed and the results are 
Summarized over the set of soundings. 

A measure of the effective cloud-cover area which has been found 
useful in previous radiational studies has been the total opaque cloud 
cover, CL, referring to the amount of thick cloud cover overhead re- 


gardless of the level. For the cloud model presented here CL may be 


expressed as 


= a k as 
CL CL, 1s CL. CL, CL. . (2-4) 


23 





Meee lennon lLAL RADIATION 


A. THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas were developed by Sasamori (1968) for flux 
emissivities in the atmosphere associated with computations for the 
radiative balance requirements of the NCAR General Circulation Model. 
Sasamori derived the empirical emissivity formulas for water vapor 
and CO. by comparison with the theoretical values built into the 
Yamamoto Radiation Chart (1952). Tne Yamamoto cnart has proved to 
be guite accurate for numerical checks of the Sasamori emissivities. 
This chart was also used in the previous Studies (A, B, C, D) as a 
systematic guide for integration of the radiative transfer formulas 
developed by Martin (1972, 1975), who adapted the Sasamori emissiv- 
ity formulas to the particular layers of interest in the gridpoint 
computations of the FNWC primitive equation model (Fig. 2). 


The essential long-wave (IR) net-flux parameters required for 


use in this study are the following: 


* 


ry 
I 


IR net flux at earth, k = 10 


10 
* 

Py =r net Elux at level k= 8 
* 

Pe = IR net flux at level k = 6 
*« 

Pi = IR net flux at level k = 4 
*& 

FP, = IR net flux at level k = 2 


In addition the IR net-flux divergence coolings to be computed at 
each gridpoint are 


Fel0 


IR net-flux divergence in the layer (8,10) 


F68 


IR net-flux divergence in the layer (6,8) 
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F46 = IR net-flux divergence in the layer (4,6) 

F24 = IR net-flux divergence in the layer (2,4). 

In the Radiation Balance Studies in the series A, B, C, D only 
F610 and F26 were computed because time restraints in the present 
FNWC Operational heating package have prevented the use of greater 
resolution in the vertical. Here, four flux divergences are computed 
in order to examine more closely the variability of the flux diver- 
gences over the layer thicknesses reduced to approximately 200 mb 
each. To compute these four flux divergences it was necessary to 


* * 
Meili ze additional formulas for F and F as developed by Martin 


8 4 
1975 )we 
In order to make IR net-flux calculations along the path of in- 
tegration, there must be a physically sound representation of the 
emissivity (coe) as a function of both water vapor and CO, absorber 
masses in layers along the sounding. For a complete discussion of 


the emissivity formulas used in the quadrature scheme, refer to 


Spaeth's Appendix A (1975). 


B. NET FLUX FORMULATIONS 

The radiative sounding as depicted in Table I was computed as 
the combination of parameters U(k,10), C(k,10) and Ty for each re- 
Gumred level, k = 10, 8, ... 1,0. Cloud parameters CL. and CL. were 


1 2 
also computed by Eq. (2-2) at each gridpoint and in general are both 
non-zero. The grid area was then considered to be composed of areal 
fractions (weights) defined in Eqs. (2-3a,b,c,d) and denoted by the 
Svnisors W(0,0), W(1,1), W(1,0), W(0,1). 


The composite net flux F (CL, CL.) at level k = 10 at each 


10 


gridpoint is then constructed by using the appropriate weight factors 
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* 
to mubeciply the reference net flux Flo computations defined for the 


four special cloud-cover cases already defined in Section II.C: 


* * * * 
E (0,0), F (1,0), F oil); F (1,1) 


10 10 10 10 


It therefore follows that 


k * * 
E (CL), CL.) = W(0,O0)F) (O70) s= W(1,0)F, (1) 


10 0 O 


== W(0,1)F) *(0,1) 5 W(1,1)P, Lene (3=1) 


0 0 


The reference net fluxes ae of Eq. (3-1) are associated with (1) 
clear skies in both layers, (2) overcast in the upper layer only, 
(3) overcast in the lower layer only and (4) overcast in both layers 
respectively. 

Spaeth (1975) has listed these reference net flux formulations 
in his Eqs. (3-6), (3-7) and (3-8). Using the definitions of W(0,0), 


aeiy@), W(O,1) and W(1,1), F..* (CL, ; CL.) can be shown to assume the 


10 


form 


* 
F (CL, CL.) SSCS 1{(B ) 


10 )-.5{€ (8,10) (B, .-B 


NOs 8 


te (EL Ney MO) E> (6710) ) (By-B,) 1} 


+ 


(1-CLy )(1-CLy ){B,-.5{(e, (6,10) 

+ Eg (4110)) (BL-B,) + (€. (4,10) 

+ € (2,10)) (B,-B,) + (€, (2,10) 

thie (CLO) WB 85 a Emme tO) 1 )*38 51) 


+ Cho ree eae (97 10) 1) ; fee) 


Here 


-/ 4 
mrs e & 
B ia O40 3 x 10 t (S—2)) 
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is the Stefan-Boltzmann blackbody flux in langlies per day. 


Further, eye WU Te 10) is the combined water-vapor and CO. emissiv- 


ke 2 


teyealtong the path from level 10 to level k. This emissivity is con- 


Sidered by Sasamori to be temperature independent for T > 210K, whereas 


=e represents the temperature dependent emissivity applicable for 


feemeclOK. [See pp. 136-137, Spaeth (1975); Sasamori (1968)]. 


* 
The formulas for Fy (CL, , CL), k = 2,4,6, and 8 have been developed 


by Martin (1975) in a manner analogous to the derivation of the weighted 


* 


Poo . The results are reproduced as the following equations: 


* 
ee Li-cr {1B,--5[€ (6,8) (BL-B ) 


8 6 


+ (E (6,8) + Ee (478)) (Be-B,) 
(Ean) a Ba ae erate) 
pacer (c co) etc yen, 8) (BSB) Sogo) meal aed | 
+ CL, [1-.5€ (6,8) ](B,-B,)+ CL 


(1-CL,) [1-.5€ (8,10)](B, ,-B ) 


1 i 8 
+ (1-CL,) (1-CL,) [1-.5€, (8,10) ] (B, )~Bg) (3-4) 
* 
nS [1-ch, ] {B.-.5le, (6,8) (B,-B,) oe OME 5) 


+ (E (4,6) + Eg (210) ) (B,-B.) + (€ (278) 
+ € ,(1,6)) (B,-B)) + € ((0,6),T,)*B,]} 

cae acupmelcl (sy be) [1-—-5(c » (6,8) 

+ S(O) LOD fein CB 18,2)" 


tec eons)]s nel crer) Wleresa 


i 
[1-,5(E, (6,8) oF Eo (Or 10) ) 13 ° (3-5) 


2) 





F, Se acbea ole (4,6) (B--B,) 
+ (€, (476) + Eo (478)) (Be-Be) + Eg (274) (BA-B.) 
+ (€ (2,4) + € (1,4)) (BL-B)) + Caen COP em =o 
+ CL, (Bye otee (2,4) (8) —B>) (3-6) 
+ (€ (2,4) + €(1,4)) (B,-B,) + ore) +B, 1s 


i (1-CL, ) (1-CL,) 0 r= pO eed CCD a Eo (4210))13 


F, = {1-cL,] iB,~-5[€, (2,4) (B,-B.) sa CN 2) 

eee) MEE aet (en (2,6) 4 _(2,8))* 

(Bg) 2s Bt 2) ee) 2 Ore) pW) ee (3-7) 
me (ilssteit Ris eine) Ea) lilo S08, (2nd) 


+ €¢ (2710) 1} + CL {B 


1 g [le (2,4) (B,-Bo) 


he (1,2) Ey ee ee, 2) cee) a2 


As was described by Spaeth (1975, Section III.B.5.) concerning the 
use of the composite case, Eq. (3-2) can be reduced to give expressions 


* 
EOD Foo for the various reference cloud-cover cases (0,0), (1,1), (1,0) 


x * 
and (0,1). The resulting schematics in the case of Fe and Fy are 


G@errecea aS the unhatched area in Figs. 3(a,b,c,d) and 4(a,b,c,d), 


* 


* 
Fe and F. can be found 


* 
respectively below. Similar graphs for Foo ; 


on pp. 41-43 of Spaeth (1975), and are not reproduced here. 


* 


A typical gridpoint listingof the IR net-flux computations, PF 


Gz=—910,0,...2), has been reproduced in Table If for gridpoint (1,1) 
based upon the radiative sounding of 16 April 1974 (see Table I). The 
printout procedure involves computation of the reference net-flux 


values 


* * * 
eos) Eo AL, 1) 


* 
me Ou, F e k 


k 
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aoe ienmometas. 9 (5—2)5,5(5-4), (3-5), (3-6) and (3-7) for level 

k= 910, 8,...2 respectively. Then Eq. (3-1) with the appropriate 
weight-factors of Eqs. (2-3a,b,c,d) has been utilized to derive the 
composite F. (CL,, CL.) Values that are listed on the bottom line of 
eeaote IT. 

TABLE II. A sample listed of IR net-flux computations, weighting 
factors and composite values, F, (CL,, CL), as computed for grid- 


Pomme (1,1) for 16 April 1974. Net flux Galues in ly min71 


Cloud Case Weight 


* * * rp” * 
F 
(CL, CL.) W(CL, , cL.) PG z F 4 i 
(O,Q) PlOe2 . 1664 £2288 . 2649 23160 ~ 3487 
laley .1497 .0943 Pluss .O915 . 1706 2338 
(9,1) 23007 HOS 12 ~1563 BO SAT Dee AD - 3202 
crt) ~4494 OS be -0586 .0607 .1706 2338 
* 
Pee composite values - 0692 lel? uot? 5 AOE SAMS 
nH 


* 


eee rROPOSPHERIC COOLING BY LAYERS; Fio COOLING 
For each mid-seasonal day listed and at each gridpoint, an IR net- 


flux computation in the format of Table II is easily converted into 


four sets of layer cooling effects. 


ae * * 
F810 = Fo = Foo 
x * 
F68 = Fe = Fo 
r * # 
F46 = Fy = Fe 
2 : ‘ 8 
F24 = FE. = Fy (3-8) 
=k 


These layer cooling rates (ly min have then been collected in meri- 
dional cross-section format for each longitude under study (see Section 


V) and by season. 
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The overall tropospheric cooling rate by IR net flux is then given 


* 


* 
simply by F, - Foo at each gridpoint for the date of the radiational 


sounding. The tropospheric cooling rates computed are then identical 


to 


* * 


- = + + + : - 
F, Foo F24 F46 F68 F810 (3-9) 


The values of F., - Flo so deduced are discussed on both a seasonal 


and a zonally-averaged basis in Section V. 
It will suffice to discuss here the zonally-averaged values of 


* 
Foo (CL), CL.) as computed by the long-wave radiational model previously 


* 


presented in this section. These results, listed simply as Flo 11391 


Table III, will be discussed as a function of seasons, latitudes and 


* 


CL (total opaque cloud cover given by Eq. (2-4)). The listings of Fao 


in Table III are essentially as extracted from computations in the 


* 


format of Table II followed by meridional-averaging of Foo across 


constant latitude lines. Finally the zonally-averaged annual values 

of both Flo and CL have been computed by arithmetic-averaging over 

the four mid-seasonal results at each five-degree increment of latitude 
meemeeOs tO GON (cf., Eq. (5-1)). 

The model-annual values of ae presented in Table III are presented 
with those derived from Budyko (1956), which in turn are listed in the 
final column of Table III. Corresponding values of total opaque cloud 
cover, CL, for the Budyko climatology were not available so that only 


a general comparison of the two annual F zonally~averaged distribu- 


10 


tions is possible. 


a2 
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* 
fable iiidepices the zonally-distributed values of Foo andeeot CL 


at the earth's surface. In a seasonal comparison of the model-computed 


* 


F values it is clearly shown that F 


* 
ss 19 (Ch) is a decreasing function 


of cloud cover. There is a clear-cut tendency in each season for a 


* 
maximum value of Foo to be located in the subtropics (latitudes 15N- 


* 


25N). Also there is evidence of a high latitude (55N-60N) minimum Foo 


associated with a concentration of maximum cloud cover CL. An outstand- 
ing variation is the transition in the Southern Hemisphere latitudes 
(20S - 10S), which has small cloud cover in local winter and compari- 
tively large cloud cover during the other three data periods. This 
cloud-cover variation corresponds in general to the ITCZ behavior in 


* 


these latitudes across the indicated seasons; so that, Foo is a maximum 


in mid-July and a relative minimum in the period January-April. 
The final entry in Table III, namely "Wt. Avg.", is the Northern 
Hemisphere mean, cosine-weighted with respect to latitude (cf., Eq. (6-1)). 


* 
The cosine-weighted F values show a minimum in Northern Hemisphere 


10 
Summer with no clear-cut differences in CL. This summer minimum is 


attributable to higher downward flux with increased summer vapor 


pressures. 


De OUTGOING IR NET FLUX TO SPACE 


1. Parameterization Formula for epee AcMosvichestkR Net biux 


Mmcnesnour previous radiat2on studies (A,.B, C, D), an approxi- 
mation to the IR net flux to space, designated as FF2, was computed 


from the radiative sounding at each gridpoint. FF2 was essentially 


* 


an extrapolation of F, to the top of the atmosphere obtained by delet- 


ing the downward IR flux due to stratospheric water-vapor and CO A 


5° 
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MoEemeEeeilse expression for the net IR flux to Space was introduced here 
after Martin (1975). The development is analogous to the quadrature 


summations 


F_ = : € dB (3-10) 


for the various reference-cloud combinations (0,0), (1,0), (0,1) and (1,1) 
of Fig. 5{(a,b,c,d) respectively. The final quadrature-formula is then 
obtained as the weighted net-flux result as was also done in Egs. (3-4), 


(3-5), (3-6), (3-7) and is listed below. 


2 [1-CL, ]{B,-.5{(e, (0,1) + 2, Ole) ae) 
+ (E072) - Eg (004) ) (B,-Bo) 


ee A072) Ei (0, 6)) (BE-B,) BE Or Sys E  (018)) (B-B 1s 


x cL, {B,--5[(e, (0,1) + Eo (0r2)) (B,-B,) (3-11) 
ae (E00, 2) fi E (074)) (B,-B,) 1h 
+ (1-CL,) (1-CL,) {{1-.5(€, (0,8) fe =, oy LO} eel 


where CLi, CL. are given by Eq. (2-2). Note that Fo Staban. (s-ll)o as 


2 
representable by the unhatched areas of Figs. 5(a,b,c,d) for the various 
reference-cloud cases and that Pe has no downward IR flux corresponding 
to the level k = 0O. 
2. Comparisons of Pie with Mid-Seasonal Satellite Climatology 
Comparison was made of computed-model values of Pin with satellite 
measurements of total long-wave flux to space (after Raschke et al., 1973) 


for the Nimbus III mid-seasonal periods most nearly comparable to that 


of the FNWC data. Pn for each mid-seasonal date was computed over the 
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four meridians considered. These Fin values were then averaged across 
the four meridians to get a mean zonal distribution of the type shown 
an Table IV. | 

Table IV shows the zonally-averaged model values compared with 
those extracted from Raschke. Raschke's results were obtained by 
averaging across the same oceanic meridians as those used in this study. 
Again, in the bottom line of each column in Table IV is listed the 
cosine-weighted mean of each set of column values for the Northern 
Hemisphere only. 

The zonally-averaged values computed by the FP model are very 
close to those reported by Raschke, especially between O-65N and in 
the Northern Hemisphericmeans. The limitations of the comparisons 
made here are obvious, when it is recalled that between latitudes 20S- 
5S and between GON-65N there are fewer than four meridional lines 
available for computing the listed zonal values in Table IV. For all 
other zonally-averaged values, four meridional lines were used in the 
averaging. 

The close comparison between model-values of a from (3-11) and those 
essentially derived from satellite climatology tend to support the 


cloud parameterization, Eq. (2-2), insofar as IR net flux is concerned. 


a 
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IV. SOLAR RADIATION 


A. COMPOSITION OF SOLAR INSOLATION 
At the top of the atmosphere (k=0) this study assumed a solar con- 
stant of 2.00 ly Phas (Joseph, 1971). Furthermore, this constant was 
assumed subject to a four percent attenuation above the tropopause due 
to ozone and oxygen. Thus the effective solar constant at level k=2 
' Se 
ime nis study is 1.92 ly min . 
To compute the effective solar insolation at the tropopause the 
following formula was used 
1 eee 
F(2) = =) ese, Cos Z (4-1) 
m 
where S is the effective solar constant at level k=2 and 
r/x = ratio of the actual earth-sun distance to the mean 
earth-sun distance, a function of the Julian date 
Cos Z= cosine of the zenith angle, a function of the 


Julian date determined by 


Cos Z = Sin @ Sin 6 + Cos ¢ Cos 6 Cosh. (A= 2) 


The symbols on the right side of Eq. (4-2) are defined as follows: 


latitude 


“a 
{I 


Omen sollaredectmnation 


25 
ll 


hour angle 
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The Smithsonian Meteorological Tables (List, 


Hm958) list the ratio of 


x/x and the solar declination, 6, for the mid-seasonal dates appli- 


cable to this study and reproduced in Table V. 


TABLE V. Values of the ratio of the earth-sun radius vector, v/r is 


and of the solar declination angle, 6, used in this study for year 


meys- 1974, 


Date r/X 
16 January 0.98372 
ie, April 2.00333 
16 July 1.01644 
16 October O299717 


6 
Bie SL Ss: 
8.48333°N 
21.50000°N 


Sa 2Z2500°S 


The value of sin ¢? was calculated using one of two different formulas, 


depending on the data-line used for the computations, in terms of the 


FNWC map coordinates (I,J) as in Eq. 


lines one may solve for I in terms of sin ¢ as in (4-3c,d): 


( 


comin 4 eee. —ae2 
97 3.752 
Line 2 Sin ¢ = ete 
9 oie (37 
ines 1,3,4 I = 32 - 22.065 
Line 2 fee e32.— 3.205 
eeietey 2) 
To St Cet Zo 
(SSP ee 88! 


=o 


ete 


4-3a,b). Conversely for these 


2(32-1) 7 
+ 2(32-1)7 


(32-1) * 
2 
+ (32-T) 


Cos 


IF + Sin 


Cos 


LF + Sin ¢ 


] 
for Lines 132 


for Line 3 


for Line 4 


(4-3a) 


(4-3b) 


(4-3c) 


(4-3d) 





Here I is the abscissa distance on the FNWC grid (Fig. 1) and varies 
aemnemadc Gescribed in Section Il. The soundings for lines 1, 2, and 
3 were all taken at OOOOGMT with the solar noon existing at the 180th 
meridian; therefore the hour angles for these three lines were 55°, 
10°, and 35°, respectively. For line 4, the soundings were taken 12 
hours earlier with solar noon at the Greenwich meridian, giving an 
Nmewenangle for line 4 of 35°. 

A very simple partition of solar insolation was utilized in this 
study after Joseph (1971). It consisted of dividing the insolation 
F(2) into two parts at level k=2. One part was considered to include 
all wavelengths » > .9 Um where absorption by water vapor and carbon 
dioxide bands are the most prevalent attenuation processes in clear air. 
This part of the solar spectrum was termed the F(A) energy and con- 
Sidered subject to water-vapor absorption but not to Rayleigh scatter- 
ing. For shorter wavelengths, A < -9 Um, absorption of the solar in- 
solation energy by water vapor was considered negligible. This part of 
the solar insolation was denoted F(S) suggestive of the fact that it 
was subject only to Rayleigh scattering attenuation in clear air. The 


two partitions are formulated after Joseph (1971) as follows: 


F (A) so49uF (2) (4-4) 


F(S) 


OO sae) (4-5) 


In this study, the introduction of two cloud decks produced cloud-reflec- 
tivity effects upon both the F(A) and F(S) solar energy insolations. 
However, in the clear areas around any gridpoint only the absorption-~ 
attenuation applies to the F(A) insolation, while only Rayleigh scatter- 


ing-attenuation applies to the F(S) insolation. 
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Be SDPLSPOSITION OF F(S) INSOLATION 
In the disposition of the F(S) insolation, Joseph (1971) determined 
maate Rayleigh scattering reflectance to space by clear skies (after 
Coulson, 1959) could be effectively approximated by least squares in 
the following form 
a(R) = .085 + .25074 [log (5~) Sec 2)] (4-6) 
O 


where Po =slOlsnzo Mon In Eq: (4-6), 1/P . =—heanewmew of the tact that 


mean sea level pressure T is close to 1000 mb. Also 
=i 
Sec z = (Cos 2) 


with Cos z given by Eq. (4-2). 
The surface albedo a(G) is another reflective parameter utilized in 
this study. Over oceanic areas the following formula for a(G) after 


Gates et al (1971), was utilized: 
gi(G) = max 1.06, .06 + .54 (.7 = Cos z)}. (7) 


As described in Section III, four combinations of reference~cloud 
cases are possible with a two-layer cloud model. The disposition of 
F(S) under each of these cases will be discussed in the remainder of 
this subsection. 

1. Clear Sky CaSe 

In the clear sky (0,0) case the F(S) insolation was subjected 
to both Rayleigh scattering reflectance a(R) and the surface reflectance 


a(G). Considering the likelihood of a succession of multiple reflections 
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between earth and atmosphere, the F(S) insolation actually penetrating 


the earth's surface after scattering is given by 


IS10(0,0) = F(S) [1-0(R)] [1+a(R)a(G) +... (a(R)a(G))™ 


fameenare | od CG) ) (4-8a) 
emeeeas, by 
m520(0,0) = F(S) [1-a(R)] [l-a(G)]/[1-a(R)o (G) ] (4-8b) 


ee = Cloudy-Sky Cases 


In the three cases in which clouds were present, F(S) insolation 
absorbed by the ground at each gridpoint was computed using the follow- 


ing equation (after Arakawa, 1972): 


menotl 1) = F(S) (1-R(1)) (1—R(2)) (1-a(G) ) 
*{1-[R(1)R(2) + R(2)a(G) + R(1)a(G) 


oa IRA aeaetenl =e (4-9) 


As indicated by the notation (1,1), denoted CL, = CL. = 1.0, Eq. (4-9) 
is the formula used in calculating F(S) insolation absorbed by the 
ground in the case where overcast clouds are present at both levels of 
Fig. 2. Also in Eq. (4-9), initial values of cloud-reflectance were 
chosen, namely R(1) = .54 for the mid-level clouds between k=4 and 6, 
and R(2) = .66 for the low-level clouds between k=8 and 9. Both cloud- 
reflectance values are as suggested by C. D. Rodgers (1967). 


For all the other cloud cases, the following changes were 


applied to Eq. (4-9). In the (1,0) case (CL, = 1.0, CL. =—"O20)- tne 
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desired earth-absorbed insolation is obtained by setting R(2) = 0.0 
in (4-9), from which it follows that 
Perot, O)y— (Ss) (1-R(1)) (1-0 (G))7/(1-R(1)a(c)]. (4-10) 


In the case (0,1), one sets R(1) = 0.0 in (4-9) so that (4-9) simpli- 


fies to 
7S10(0,1) = F(S) (1-R(2)) (1-a(G))/[1-R(2)a(G)]. (4-11) 


Note that with a cloud overcast present, the Rayleigh clear-sky scatter- 
ing @(R) does not appear in Eqs. (4-9), (4-10) or (4-11), but is included 
empirically in the cloud reflectances R(1) and R(2). 

sm Composite F(S) Insolatien at Earth 

Equations (4-8), (4-9), (4-10) and (4-11) were utilized in the 

computation of the cloud-weighted F(S) insolation penetrating the 
earth's surface considering the areal-weights of the cloud combinations 
denoted by (0,0), (1,1), (1,0) and (0,1) about a gridpoint. The re- 
sultant F(S) insolation penetrating the earth's surface denoted by 


ISl10 is therefore expressible as 


IS10(CL,, CL,) = 1810(0,0) W(0,0) 
TEGO, 1) ae) 
STO Onna o) 


+IS10(0,1) W(0,1) . CID) 


Here the weighting factors W(0,0), W(1,1), W(1,0) and W(0,1) are com- 
puted in Eqs. (2-3a,b,c,d) respectively. Note finally that the part 
of F(S) insolation reflected to space is found by subtracting IS10 


(CL, , cL.) EBON GE CS)". 
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Table VI lists the F(S)-disposition resulting from a particular radia- 
Eivewsouncding at gridpoint (1,1) on 16 April 1974. The individual 
computations of IS10 as they apply for the possible overcast-clear 


layer cases are made under the heading "IS10." The difference 
Rubs = F(S) = JSl0 (4-13) 


in each case represents F(S)-insolation reflected to space while 


IS10 
STRAN = —T—a7(G) (4-14) 


has been computed as that portion of the F(S)-insolation incident at 
the sea surface just prior to transmission by the surface. Note that 
no absorption in air has been included in the computations of Table VI, 
and that the only absorption permitted is that implicit in IS10. 
Finally at the bottom of each column, e.g., IS10, the composite value 
has been computed by means of the weighting scheme of Eq. (4-12). 

TABLE VI. A sample listing of values of F(S) insolation (ly ain} 


Comeuued at gridpoint (1,1) for 16 April 1974 using Eqs. (4-6,..., 
4-14). 


Cloud-case Weight 

(CL, , CL.) W(CL CL.) £510 REEFS STRAN 
(OF 0) -0993 s43105 e) solo so2 lo 
GIL 10) ~ 1484 SAE, SS ly: (SOAS 
(0,1) 6 SCANS. Seal ono e2327 
C18) ASO 1340 -4716 - 1696 
Composite-F(S) values -2014 ~4041 ~ 244] 
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C. DISPOSITION OF F(A) INSOLATION 
The fractional portion of the solar insolation subject to absorption 
by atmospheric water-vapor and carbon dioxide 1s covered in the follow- 
ing subsections. 
1. Clear-sky Case (0,0) 
The Manabe-Moller absorptivity function provided the necessary 
absorptivity values for the key layers in this case. The form of this 


absorptivity function is 


a(2,k) = .271(U(2,k) Sec 2]? . (4-15) 


PememaetS the absorptivity applied to the pressure-scaled water vapor 
mass between levels 2 and k (Fig. 2) along the zenith slant-path angle 
z. The resultant absorbed insolational energy in the particular layer 


(2,4) is then given by the Manabe-Moller relation 


A24 = 0.271F (A) (U(2,4) Sec Sie es ; (4-16) 


In the same manner A26, A28 and A210 are found. Then the absorbed in- 


solation in the layers (4,6), (6,8) and (8,10) are computed by 


A46 = A26 - A24 (4-17) 
A68 = A28 - A26 (4-18) 
A810= A210- A28 . (4-19) 


Water-vapor mass above level 2 was assumed negligible in the F(A) dis- 
position of the solar insolation. 

By subtracting A210 from F(A), the direct transmission of F(A) 
insolation impinging at the earth's surface was determined. The trans- 


mission of F(A) insolation is then further reduced by the transmissivity 





(l=a(G)) after surface-reflectance,which leads to the earth-absorbed 


mas@dlation 


Zoe 


TA10(0,0) = F(A) {1-.271[(U(2,10)Sec 2z] Glo. (G) )- (4220) 


The transmitted energy impinging upon the earth just prior to absorp- 


aren 1S 


TRANA(O,0) = IA10(0,0)/[1l-a(G)] . (4-21) 


In the remainder of this subsection, which discusses the cloudy 
layer cases, representative cloud reflectivities and cloud absorptivi- 
Mmreeewere initially adopted, after C. D. Rodgers (1967), for the two 
possible cloud layers. These initial cloud reflectivities were RA(l1) = 
-46 and RA(2) = .50 while the cloud absorptivities were A(1) = .20 and 
A(2) = .30. Note that the reflectivities for the F(A) wavelengths are 
somewhat smaller than those adopted for the F(S) wavelengths. The pro- 
cedure of considering the cloud conditions to be overcast whenever they 
appear and then applying the appropriate weighting factors in the com- 
posite summation will again be followed as in Sec. IV.B. 

To simplify the discussion for the cloud-covered cases, the (1,1) 
case will be presented first as it contains representative type equa- 
tions for the remaining two cases, (1,0) and (0,1). 

2. Overcast in Both High- and Low-cloud Layers 

The following set of formulas illustrate the model disposition 
of incoming solar insolation (F(A)) from level k=2 to the earth's sur- 
face and permits determination of the amount of insolation absorbed by 
the atmospheric layers and by the earth's surface. The dashed separa- 


tion lines are introduced to subdivide the absorption and reflection 
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physics of the model into subsections which permit the analysis to 
proceed more or less within successive 200 mb layers. The equations 
relate to the parameters in Fig. 6 where the insolations, A24, Ad6, 
A68, A89 and A910, etc., represent the contributions to the insola- 
tion absorbed in the layers involved. Symbols F2, F4, F6, F8 and F9, 
etc., depict the streams of insolation passing through the indicated 
level. A vertical arrow implies the direction of insolation passage, 
i.e., ¥ denotes downward insolation, t upward-reflected insolation, 
and ++ downward-reflected insolation. Terms involving the symbols 
"TD", as expressed by the functions of Eq. (4-22), indicate the Manabe- 
Moller transmissivities for diffuse insolation in the layer beneath an 
Swstang Cloud. In the latter situation, the term Sec Z in Eq. (4-15) 
is effectively replaced by the mean slant-path, eeeua = 5/3 (Katayama, 


ico), €.9., 


3 


TD68 = 1-.271[U(6,8) 5/3]°>° 


STANT ONS aie “ (4-22) 


EDI LO 


In accordance with the above definitions and as described by Fig. 6, 


the formulas for the (1,1) case are listed: 


F4+ = F(A) (1-.271[U(2,4) Sec ies) 

F4t = F4+(RA(1)) 

F2t¢ = F4t (1-.271[U(2,4) Sec 2 eas 

A24 = F(A) - F4t) + F4t - F2t (4-23a) 
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A46 = 


F4y (A(1)) 


me wm em i ee SS SO SP SS SS ll ES ll OPO. Ell el 


F6Y = 


F8Y = 


F8¢t = 


F6t = 


F6y+= 


F8Y¥= 


A68 = 


F4y - F4t - A46 
F6¥ (TD68) 

F8¥ (RA (2) ) 

F8+ (TD68) 

F6t (RA(1)) 

F6+¥ (TD68) 


F6Y - F8Y + F8t -— F6t + F6YY - F8YY¥ 


A89 = 


FOY = 


F1l0¥= 


F10t= 


pov = 


Eoyy 


F1l0VY¥= 


A910 


A810 


Note 


(F8Y + F8¥¥) (A(2)) 

F8¥ - F8t + F8+¥¥ - A89 
F9} (TD910) 

F1OY (a (G)) 

EL OAR S10) 

F9*+ (RA(2) ) 


Bo Do LO) 


FOY - F1lOY + F1Ot - FO*V + FOVY - F1OYY 


A89 + A910 


(4-23d) 


that the effect of multiple reflections between clouds or 


between the earth's surface and the lower cloud has been incorporated 


to include only the effect of two reflections, with the lowermost re- 


flecting surface absorbing the remaining impinging insolation. 


Computa-— 


tions indicated that the insolation remaining after two reflections was 
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F(A) 

. Bo 

A24t A24y 
4 
2 RA (1) 

F6oty+ F6t F64: 

A68Yy Aest A68y 

8 

A910¥+ A910+t A910 
7 EEOt F10¥ a (G) 

Jo Ge ON TV VI a a ee eA 
IA1LO 


Figure 6. Schematic representation of F(A) insolation disposition in 
the case of two overcast layers, 
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too small to warrant the consideration of further reflections. Also 
insolation reflected upward from a lower interface (cloud or ground) 
to the base of an upper cloud deck has not been subjected to upper 
cloud absorption. This tends to reduce very slightly the secondary 
cloud-absorption. 

From Eq. (4-23d), the impinging F(A) insolation at the earth's 


surface can be expressed 


TRANA (1,1) = F1OY + F1OYVY (4-24) 


The F(A) insolation which 1S actually absorbed by the earth's surface 


(see Fig. 6) may be written as 


IA10O(1,1) = F10¥(1-a(G)) + F1OVY (4-25) 


eS. Drsposition Cet Caelnsolatwon with an Upper Overcast Only 


With a single cloud layer present only at the upper level (Fig. 
7) the equations depicting the model-disposition of incoming insolation 


becomes a simplified subset of the previous case: 


F4yY = F(A) (1-.271[(U(2,4) Sec Bie 

F4t+ = F4Y (RA(1)) 

Pete—erot(l—.271 (002.4) Séc ie 

A24 = F(A) —- F4¥ + F4t —- F24 (4—26a) 
Ado = F441) (4—26b) 


awh 





F(A) 
F2+t 
2 
Feat A241 
F4t RA (1) F4y 
4 
PoC) A46 
6 er ree 
RA (1) 
Footy 6t F6y+ 
U | 
A68ty+ A6st A68y 
8 
F8t¥ PS8t F8Y 
A810+¥ A810+ A810 
F10t F10¥ a (G) 
10 
Vee 7 7 SLAP SS SSS 


IA10 


Figure 7. Schematic representation at F(A) insolation disposition 
with an upper overcast layer only. 
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F6Y = F4) - F4t - A46 
F8¥ = F6+(TDE8) 

FlOY = F6¥(TD61LO) 

FlOt = F1O¥(a(G)) 

F8t = F10t (TD810) 

F6t = F10¥ (TD610) 
F6¥¥ = FO6t(RA(1)) 


F8\) = F6+Y¥ (TD68) 


FlLOYVY = F8Y)(TD810) 
A68 = F6t - F8Y + F8St - F6t + F6YY - F8YY (4-26c) 
A810 = F&Y - FlOY + FlOt - FSt + F8YY - FI1OYVY (4-26d) 


The variables used above are defined in a similar manner to those in the 
(1,1) case. The impinging insolation at the earth's surface may be for- 


mulated as follows 


TRANA(1,0) = FlOY + FlOVY (4-27) 


while the F(A) insolation absorbed by the earth in this case is given by 


TA1O(1,0) = F10¥(1-a(G)) + FlOVY . (4-28) 


aS 





4. Disposition of F(A) Insolation with a Low Overcast Only 


With an overcast lower cloud layer the model-disposition symbols 


are aS depicted 


F44 = 
F6+ = 
F8¥ = 
F8t = 
Fot = 
F4t =: 
F2¢ = 


A24 = 


in Fig. 8 and are physically related as follows: 


POMRGIen o7i(2.4) Sec 21°”) 
Pyles (2,6) Sec mo 
Biba 27)1U(2,68) Sec ns 
F8¥ (RA(2) ) 

F8t(1-.271[U(6,8) Sec ayes 
Fat (1-.271[U(4,8) Sec 2]°>">) 


F8t(1-.271[U(2,8) Sec z]°°°7) 


F(A) - F4+ + F4t - F2+ 


= we wee is ee ee i ss i eee ee i lee 


A89 = 


no, = 


FlO¥ = 


Fl1Ot = 


FOF = 


FOV} = 


F1l0++ = 


A910 = 


A810 = 


F8Y (A(2)) 

F8+ - F8t - A89 

F9+ (TD910) 

F104 (a (G) ) 

F10%+ (TD910) 

FOt (RA(2)) 

F9+¥ (TD910) 

FOY - F1lOY + F1lOt - FOt + FOYY - F1OYVY 


A89 + A910 
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F(A) 
F2+ 
2 
A24+ A24y 
FAt FP4y 
A46+ A464 
6 
F6t F64 
A68t A68Y 
F8t RA (2) F8Y 
8 | 
A(2) A89 
FOV) , F9t RA (2) FY 
A910+4 A910+ A9104 


rt | | | F10¥ 
ee el a i ff ff f/ /// 


IA1O 


Figure 8. Schematic representation of F(A) insolation disposition 
with a lower overcast layer only. 
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Again the variables are defined as stated before in the (1,1) 


Case. 


MRANA (OQ, 1) 


while that portion which is absorbed by the earth (Fig. 8) is 


IA10(0,1) 


ineeqs. (4-25), 


= FIO)’ + FlOYVY 


= F10¥(1-a(G)) + F1lOYY 


> 


piee wise nesinetcent tlux 1S defined at the earth's surface as 


(4-30) 


(4-31) 


(4-28), and (4-31), the quantity F104+ is small enough 


in each case so that no further reflections from the earth were con- 


Sidered. 


5- Composite F(A) Layer-Absorptions and Surface-Absorption 


Insolation 


As has been previously discussed, the standard grid-area weight- 


ing scheme of this study was applied to obtain composite values of the 


absorbed F(A)-insolation in key layers and also within the earth's sur- 


face. 
combination absorption 


results: 


A24(CL,,CL,) = 


A46 (CL, ,CL.) = 


A46(CL, ,CL.) = 


The weighting factors applied to the corresponding overcast- 


quantities provided the following composite 


A24 (0,0) 


A24 (0,1) 


A46 (0,0) 


A46 (0,1) 


A68 (0,0) 


ASS (Om) 


W (0,0) 


W(0,1) 


W(0,0) 


W (0,1) 


W(0,0) 


W (0,1) 
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A24 (1,0) 


A24(1,1) 


A46 (1,0) 


A46(1,1) 


A68 (1,0) 


A68(1,1) 


W(1,0) 


W(1,1) 


W(1,0) 


W(1,1) 


wil? 0) 


Walt) 


(4-32) 


(4-33) 


(4-34) 





t! 


A810(CL, ,CL,) Polomenojewto,O) + ASIO(1,0) W(1,0) 


A810(0,1) W(0O,1) + A810)1,1) W(1,1) (4~35) 


+ 


I 


TALO(CL, ,CL.) PeOrO,Ojmw (0,0) + TAIO(1,0) W(1,0) 


= 


mle Om (O,)) + TALO(1,1) W{1,l) . (4-36) 


[Mire weighting factors W(0,0),...,(W1l,1) were listed in Eqs. (2-3a,b,c,d), 
and A24(0,0), A46(0,0), A68(0,0), A810(0,0) and IA10(0,0) are given in 
each clear sky case (0,0) about each gridpoint by Eqs. (4-15), (4-16) , (4-17) 
(4-18), (4-19) and (4-20) respectively. 

The results for the absorption in layers (2,4), (4,6), (6,8), 
(8,10) and at the surface, level k = 10, are shown in the following table: 
TABLE VII. A sample listed of F(A) insolation values (ly min ~) 
eempueced at gridpoint (1,1) for the 16 April case. 


Cloud-Case Weight 


CL), CL, W(CL,,CL,) A24 A46 A68 A810 IAlO0 =REFA TRANA 
(0,0) .0993 .0454 .0356 .0386 .0418 .1348 .0285 .1634 
(10) .1484 Roce smnOSS mE ODCONNeO156.9 20482 .1137 .0581 
(01,1) .3016 .0464 .0381 .0728 .0771 .0239 .0664 .0286 
(le )) aS O7 ROCs 3 O50 mmr O432..020555.0104 21233 .0126 


Composite-F(A) values TOs oO menOCe 5 mmnO204) 2. 0425.55.0324 .0953 .0391 


In the computational scheme indicated by the entries of Table 
VII, the reflected F(A) insolation to space has been depicted by the 


symbol REFA, and its values follow from 


REFA = F(A) —- A24 - A46 — A6S —- AS1O —- IALO (4-37) 


Si) 





whereas the TRANA dispositions are given by Egs. (4-21), (4-24), (4-27) 
and (4-30) or by its weighted-mean value in the case of TRANA-composite. 
6. Absorptivity (ABA) by Layers 
Here the (fractional) absorptivity as well as the actual insola- 
tion values absorbed in the layers are considered. In the computation 
of absorptivity, which is fractional absorption, the total undepleted 
insolation at the top (k=0) is used as a base. The following equation 


@easeutilized in this calculation: 
-—2 
FADJ = 2.00(x/r_) CosaZ (4-38) 


The absorptivity of the troposphere ABA was computed from the ratio of 
the insolation absorbed in the troposphere to the insolation incident 


at the top of the atmosphere rather than at k=2: 


_ A24 + A46 + A688 + A810 (4-39) 
oe FADJ 


D. ALBEDO (ALB) OF THE EARTH-TROPOSPHERE SYSTEM 

In considering the planetary albedo, the reflected insolations of 
the earth-troposphere systems in both the F(A) and F(S) insolational 
regions must be recalled by the program. Thus REF is computed at each 
gridpoint as the sum of the reflected insolation energy in F(A), denoted 
REFA in Eq. (4-37) and the reflected part of F(S) denoted REFS in Eq. 


(4-13): 


REF = REFS + REFA (4-40) 


Finally the planetary albedo is related to FADJ through 


ALB =.——— , (4-41) 
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a 


E. COMPOSITE ABSORPTIVITY (ABG) BY THE EARTH-SURFACE; 
COMPOSITE ATMOSPHERIC TRANSMISSIVITY (ATRAN) 


1. Absorptivity (ABG) of Earth 
By summing the weighted values of F(S) and F(A) portions of 
the incoming insolation entering the earth, the total insolation ab- 
sorbed at the earth's surface was computed. This quantity when 
divided by the extraterrestrial insolation gave the fractional absorp- 


tivity (ABG) of the earth's surface. The equation for ABG was 


TRO + ES1 0 
a FADJ oa) 


where IA10, IS10, and FADJ were defined previously by Eqs. (4-36), 
(4-12) and (4-38) respectively. 


Peeeueansmissivity (ATRAN) Of the Troposphere 


Also computed was the total insolational energy TRAN, incident 


Gemeieweadrtn's Surface just before absorption by the surface. This 


calculation is given by 


TRAN = TRANA + STRAN . (4-43) 


Here STRAN = [IS10/(1-a(G)] was previously defined in Eq. (4-14) and 

a(G) was given in Eq. (4-7). TRANA has also been defined as the weighted 
Value of TRANA (0,0), TRANA (1,1), TRANA (1,0) and TRANA (0,1) given by 
Eqs. (4-21), (4-24), (4-27) and (4-30). Note also in justification of 
STRAN that the four cases for IS10 of (4-8), (4-9), (4-10), and (4-11) 
each have the common factor (1-a(G)) in the numerator and therefore each 
transmitted F(S) insolation component available at the earth just before 


absorption needs only be divided by (1-a(G)). TRAN may thus be viewed 


femerwemeOtameinsoladtaonal energy incident at a pyyheliometer located at 


ao 





earth. The (fractional) transmissivity of the troposphere (ATRAN) is 


then computed from 


ATRAN = TRAN/FADJ . (4-44) 


Note finally that the major dispositions of the total insolation at 
the indicated map times have now been identified by the fractional 
values, ALB, ABA or ABG, and ATRAN, representing the reflectivity 
(albedo), absorptivity of air or earth, or atmospheric transmissivity 


as the case may be. 


F. ALBEDO TUNING BY COMPARISONS WITH SATELLITE CLIMATOLOGY 
1. General Remarks Concerning a Need for Tuning Albedos 

In the four previous radiational studies (A, B, C, D) compari- 
sons were made between the solar-insolation albedo-model (ALBMOD) com- 
putations and the satellite-climatology albedo (ALBRAS) of Raschke 
et al. (1973). The results of these comparisons indicated excessively 
high values of ALBMOD especially in the tropical and subtropical oceanic 
areas. The question was raised by the previous investigators as to 
Whether vertically-structured convective cloud elements in the tropics 
would have as high a reflective capability as attributed to the large- 
scale cloud masses existing primarily in horizontal layers as indicated 
by the parameterization of Eq. (2-2). It should be recalled that the 
initial cloud-reflectances used were as suggested by C. D. Rodgers 
mive7jeand are listed in Secs. IV.B.2. and IV.C.1. 

For the purpose of tuning the model solar dispositions, and 
in particular ALBMOD with ALBRAS, the initial values of Rodgers' cloud- 


reflectances were simply adjusted by a multiplicative factor (f) which 
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turns out by the least-squares fitting to be smaller than unity. This 
device allows a greater fraction of solar radiation to penetrate down- 
ward through the cloud-types in both the tropical and subtropical areas, 
and provides better agreement between ALBMOD and ALBRAS. 
2. Method of Tuning ALBMOD to ALBRAS 

The model-tuning process involved several steps. First, the 
new variable f was defined and allowed to range in value from 1.2 to 
0.25 by 0.05 intervals. The f-values were then multiplied by the 
initially used cloud-reflectance values (after Rodgers, 1967) giving 
twenty sets of new cloud-reflectivity values. For each set of new 
values, ALBMOD and subsequently the difference A(ALB) between ALBMOD 


and ALBRAS; l.e., 


A(ALB) = ALBMOD - ALBRAS (4-45) 


was recomputed at each gridpoint for the four midseasonal dates. 

Because cloud-types differ in structure over tropical and extra- 
tropical oceanic areas the resulting A(ALB) data, from (4-45), were 
divided into two populations denoted TR and HL. The population TR in- 
cluded the set of gridpoints located south of and including 25N, while 
HL included the set of gridpoints located north of 25N over all four 
meridians. 

Only those gridpoints where the value of at least one of CL 


1 


and CL. was greater than 0.1 were included in those populations. The 
gridpoints so excluded were considered climatologically unrepresenta- 


tive for the purpose of cloud-tuning the albedo. The A(ALB) data was 


additionally divided by season giving a total of eight separate cases 
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Partner defined in sequence as WIN-TR, WIN-HL, SPR-TR, SPR-HL, SUM-TR, 
SUM-HL, AUT-TR and AUT-HL. 

The Bimedical set of programs (Dixon, 1973) was utilized to 
compute the root-mean-square (RMS) of A(ALB) values corresponding to 
a wide range of f-values for each of the seasons and regions. Values 
of £ leading to the minimum RMS of A(ALB) were isolated and appear in 
Table VIII, along with the RMS difference of A(ALB). 

TABLE VIII. Values of f and corresponding minimum RMS 
values of A(ALB). 


Minimum RMS 


f of A(ALB) 
WIN-TR 40 7069 
WIN-HL 12119 2073 
SP Roun 35 OS 
Se Re mia als)S) = O}e 
SUM-TR - 60 O27 
SUM-HL a7 0 209 
AULT-TR - 40 SOD 
AUT-HL a(S 10) Oe 


As can be seen from Table VIII, the listed values of f in 
Winter indicates the largest difference between the tropical and extra- 
tropical cloud-types. The evenly layered cloud formations of the 
northern latitudes in mid-January seems more closely to approximate 
those considered by C. D. Rodgers (1967). On the other hand, the 
smallest difference in f-values occurs between SUM-TR and SUM-HL 


indicating a less significant difference in cloud types over the 
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complete latitude range during mid-July. As would be expected the 
Peeing and Autumn f-ranges fall in between the Winter and Summer cases. 
A full listing of zonally-averaged ALBMOD values, after tuning at all 
gridpoints, and of zonally-averaged ALBRAS values, generated from 
identical gridpoints, can be found in Table IX. 

Upon examination of Table IX it iS apparent that the tuned 
model-albedo values compare closely with the satellite albedos derived 
from Raschke. Table IX includes the annual averages and the Northern 
Henisphere cosine-weighted means as well as the mid-seasonal averages. 

Further inspection of the weighted means shows that ALBMOD is 
somewhat smaller than ALBRAS for each of the four seasons and annual 
average. However most of this underestimate seems to occur in the 
HL-population where the result of minimizing A(ALB) led to values of f 
substantially less than unity (Table VIII). This result was unexpected 
and is probably due to a mismatch between the model clouds (Eq. (2-2)) 
and those averaged in the satellite climatology particularly in the 
higher latitudes. 

The conclusion is that reflectance-tuning to satellite albedo 
appears valid in the tropics and subtropics where cloud-formations 
On a given day tend to be persistent. However tuning by this method in 
high latitudes is less conclusive because interdiurnal cloud variability 
1s much greater. 

In any case, all further computations of solar dispositions in 
this study will utilize the appropriate cloud-reflectance tuning factors 


as listed in Table VIII, depending upon season and latitude. 
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V. MERIDIONAL CROSS-SECTIONAL DEPICTION 
OF THE RADIATIVE~BALANCE COMPUTATIONS 
fe GENERAL 
The general design of this section is to utilize all of the compu- 
tational concepts discussed in Sections III and IV in the computations 
for a single time-step in the radiative heating package developed for 
use in the FNWC prediction model. After testing ALBMOD with the corres- 
ponding values of Raschke et al. (1973), it was decided that only the 
computations by the appropriate cloud-reflectance tuning factors (Table 
VIII) would be displayed in the meridional cross-sections (Figs. 10... 
17) and in the comparisons of the four~layer versus the two-layer cool- 
ing rates summarized in Table X. 
The radiative calculations were performed at each gridpoint of the 
four meridians for each of the four mid-seasonal days considered; how- 
ever, only the winter and summer cases, as illustrating the two extreme 


Situations, are presented in this section. 


B. GEOGRAPHICAL REPRESENTATION OF THE RADIATIVE~BALANCE DISTRIBUTION 
The FNWC gridpoint processed analyses for OOOOGMT were used at the 
three Pacific cross-sections, while that for 1200GMT were used for the 
Single Atlantic meridian. This was done so that the set of gridpoints 
was considered to be subject to the actual radiative-transfer calcula- 
tions involved for these specific times. Figure 9 depicts in symbolic 
language the key to the computational entries in Figs. 10, 1l,...16, 
17. This symbolic list presents the computations made at each radia- 


tive sounding gridpoint (I,J)having data in the form of Table I. The 
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computations proceed from the tropopause (approximately level k=2) to 
the ocean surface. For purposes of climatological data comparison, 
pagel, 11,...16, lL? were developed by interpolating gridpoint results 
to integral multiples of 5-degree latitudinal increments. The inter- 
polation routine to this gridpoint spacing made use of the Lagrangian 


@ieic interpolation scheme (after Spaeth, 1975). 


iba S ) (1-4) (in 


Q(T) = Q) ray (r-3) (7-4) * 22 (2-1) (223) (2-4) 


(5-1) 
(ie) (t=2) (1-4) it (isl) i 2a Ga), 


+ i a a a 
See Se )(5=4) 4) (4211) (4-2) (4-3) 
Finally for ease in reconciling the magnitudes of all radiative-trans- 
fer rates, the time-dependent solar disposition rates have been 


averaged to 24-hourly rates. 


C. EXPLANATION OF SYMBOLIC TERMS 
1. Cross-Section at Level k=2 
The discussion of all insolation parameters discussed previous- 
ly in Section IV dealt with the specific time of day that corresponded 
to the hour angle h for the instantaneous time t under consideration. 


The incident solar insolation dealt with is then 


F(2) = S(=)"* cos z . (se2) 
m 
In order to avoid reference to specific map times t, the instantaneous 
solar hour-angles were h = 35, 10, 55, and 35 degrees, respectively 
for cross-sections 1,2,3 and 4 as depicted in Figs. 10, 11, 12 and 13 


Ponewimteer, and in Figs. 14, 15, 16, and 17 for summer. 





OAVE represents the 24-hour average of F(2) and appears as the 
first input symbol in Fig. 9. Its value is considered to be more re- 


presentative climatologically for the data day under consideration than 














ive) 
OAVE iS derived by the formula 
AVE = F(2) OSs ie) 
2 Cos Z 
where 
Cos z = [H SindSind + CosdCosdéSinH] /7 (5-4) 
H = ArcCos [-TandTané] (5-5) 


Here 6 is the appropriate solar declination angle as listed in Table V, 
and H is the appropriate hour angle at local sunset at latitude $6. The 


value of H also depends upon which mid-SseasSonal date is being considered. 





Cos z in Eq. (5-4) is equal to the 24-hour average coSine of the zenith 
angle, Eq. (4-2). The 24-hour time averaging period for QAVE gives heat- 
ing results consistent in magnitude with the terrestrial flux divergences, 
which change only Slightly with the time of day. The conversion to ex- 
pected daily averaged solar disposition quantities is compatible with 

the determination of a hemispheric radiative balance for the given mid- 
Seasonal dates (Sec. IV.A.). 


Other parameters needed for level k=2 are 








= COse2Z ‘a 
QREF = REF(t) ((ooy) (5-6) 





where 


REF(t) = F(2) - A24 - A46 —- A68 —- A810 


- (IA1O + IS10) (5-7) 
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REF (t) is the instantaneous solar reflected insolation at a gridpoint 
and QOREF is its 24-hour average, assuming that the instantaneous 
planetary albedo remains constant for the 24-hour period. This assump- 
tion requires that the cloud amounts computed at the indicated synoptic 
times are representative of the entire day. 
The same principle will be used with regard to all other solar 
: parameters in the conversion from time-dependent values at solar time t 
to 24-hour averaged values. Superior bars ( ) are not used eee sym- 
: bolism for the averaged values shown in the cross-sections key, Fig. 9, 
but are implied by the use of QAVE, etc. in the solar-disposition terms. 


The 24-hour average tropopause balance, BALT, iS computed from 
* 
BALT = QAVE - (QREF + F. ) (5-8) 


for the level k=2 at the indicated latitude. Net terrestrial fluxes, 


* 


Such as F., , were considered to be constant throughout the 24-hour period, 


a valid assumption if the cloud cover remains quasi-constant for the 


| 
period. 
2. Cross-Sections in Layers (2,4),(4,6),(6,8) and (8,10) 
The following definitions apply for the four layers identified 
in the present section heading and are further identified in Fig. 9. 
All of the heat transfers shown in these layers are assumed to be of 


radiative character only. The daily-averaged radiative heating (cooling) 


| rate in layer (2,4) is given by 


: BAL 24 = 924 - F24 (5-9) 
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Figure 10(a). 
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125W Longitudinal cross-section, tropical section. 





computed from data for 16 January 1974. 
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125W Longitudinal cross-section, higher latitude section. 


Figure 10(b)*. 


Values computed from data for 16 January 1974. 
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Figure ll(a). 
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Refer to Fig. 


170W longitudinal cross-section, higher latitude section. 


Values computed from data for 16 January nova. 


Figure 11(b). 
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Figure 12(a). 
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145E Longitudinal cross-section, tropical section. 


computed from data for 16 January 1974. 
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Refer to Figwm@e for key. 


145E Longitudinal cross-section, higher latitude section. 


Figure 12(b). 


Values computed from data for 16 January 1974. 
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35W Longitudinal cross-section, tropical section. 


computed from data for 16 January 1974. 
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Hague 13 (b). 


Refer to Fig. 9 for key. 


35W Longitudinal cross~section, higher latitude section. 


Values computed from data for 16 January 1974. 
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Figure 14(a). 
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125W Longitudinal cross-section, tropical section. 


computed from data for 16 July 1974. 
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Figure 14(b). 


Rote lto Pig ..9 for key. 


125W Longitudinal cross-section, higher latitude section. 


Values computed from data for 16 July 1974. 
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Figure 15(a). 
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Ronee toubag, 9 for key. 


170W longitudinal cross-section, higher latitude section. 


Values computed from data for 16 July 1974. 
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Figure 16(a). 
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145E Longitudinal cross-section, tropical section. 


computed from data for 16 July 1974. 
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145E Longitudinal cross-section, higher latitude section. 


Values computed from data for 16 July 1974. 


Refer to Fig. 9 for key. 


Froure 16(b). 
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Figure 17(a). 


Values 


Refer to Fig. 9 for key. 


35W Longitudinal cross-section, tropical section. 


computed from data for 16 July 1974. 
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Reter SeOe Pad 9 tor key. 


, higher latitude section. 


35W Longitudinal cross-section 


Values computed from data for 16 July 1974 








where 


Q24 = daily-averaged solar absorption in layer (2,4) 


and is defined relative to A24(t) by a cosine transformation similar 


Mer bq. (5-3), and 


Hae = terrestrial cooling rate in (2,4). 


Similarly the daily-averaged radiative heating (cooling) rates in 


layers (4,6), (6,8) and (810) are given by 


BAL46 = Q46 - F46 (5-10) 
BAL68 = Q68 - F68 (5-11) 
BAL810= Q810- F810 (S12) 


respectively. 
3. Cross-Section at Air-Sea Interface (k=10) 
The radiative balance at the earth's surface (BALB) is as 


defined in the following equation: 


BALB = QABG - FP : (5-13) 


OABG is the 24-hour average solar insolation absorbed by the surface 


as follows: 





OABG = QABG(t) (Cos z/Cos z) . (5-14) 
D. MERIDIONAL CROSS-SECTIONS OF THE VERTICAL RADIATION BALANCE 


Eros. 10, 11,...i6, 17 as previously Gxplagmed represent the 


Single time step of heating computations for each of the four meridians 
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used in this study. The eight figures have been divided into (a) tropi- 
cal results and (b) mid-to-high-latitude results for the mid-winter and 
mid-summer cases respectively. While the results depicted in these 
cross-sections are exhibited as representing daily-averaged values, they 
are actually based upon radiation-computations at the specific map times 
of OOOOGMT and 1200GMT on 16 January and 16 July 1974. Therefore, for 
these results to be meaningful as a stepwise part of the heat package 
subroutine of FNWC, the solar radiative absorption and reflectance terms 


mena have to be recoverable as a function of GMT, i.e., 





F(2,t) = QAVE* (Cos z/Cos 2z) (5-15) 
REF(t) = QREF* (Cos z/Cos 2) (5-16) 


etc. Thus solar disposition terms may then be utilized in connection 
with the one-hour stepwise application of the thermodynamic equation of 
the set of primitive equations used in the FNWC prediction process, 


assuming the "2/3-CL" parameterization of Eq. (2-2). 


E. COMPARISON OF THE FOUR-LAYER WITH THE TWO-LAYER FLUX DIVERGENCE MODEL 

In previous studies (A, B, C, D) which considered only the two-layer 
flux-divergence model [layers (2,6) and (6,10)] it was assumed that the 
radiative-cooling rate in layer (2,6), BAL26, was uniformly distributed 
over layers (2,4) and (4,6). A Similar remark is applicable to BAL610 
relative to the sublayers (6,8) and (8,10). In an effort to determine 
the vertical resolution of the atmospheric flux-divergences separately, 
the four values of BAL24, BAL46, BAL68 and BAL810 have been averaged 


Over all gridpoints for the four computational dates. 
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The computed values of BAL24 and BAL68 were then compared with one- 
half of BAL26 and BALG1O respectively to determine the relative per- 
centage of heating (cooling) in these layers, as depicted in Table X 
for the mid-seasonal winter and summer dates. It is noted that the 
values of BAL46 and BAL810 are identical in magnitude but of opposite 
algebraic sign to BAL24 and BAL68 respectively. 

The results of these four-layer computations (Table X) indicates 
that in layer (2,4) there is approximately thirty percent less cooling 
than would have been deduced in the simpler two-layer model. This 
implies less radiationally induced instability in the upper troposphere 
than that originally estimated. In the layer (4,6) there is compensat- 
ing increased radiative de-stabilization. This radiative difference 
in the layer (4,6) appears to be primarily a result of the placement of 
clouds in both models. For instance, the upper clouds have tops at 
k=4, and IR net flux acts as the chief source of increased de-stabili- 
Zation in layer(4,6) (and of decreased de-stabilization in layer (2,4)). 
Likewise the lower clouds have tops at k=8; so, it follows that IR net 
flux acts as the chief source of increased de-stabilization in layer 
(2,10), and finally, of reduced de-stabilization in layer (6,8). 

The tuning of solar reflectances has not introduced a significant 
radiative influence on the results of the de-stabilization differences 
Just outlined. Rather it was the placement of clouds relative to 
layer boundaries that affected primarily the IR net-flux calculations. 
The results summarized in Table X are then to be considered valid if 


the cloud-parameterization model approaches climatological reality. 
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Vi- THE ZONAL DESTRISUTION OF RADEATIONAL 
BALANCE TERMS OF THE OCEAN-ATMOSPHERE SYSTEM 

A. GENERAL INTRODUCTION: ZONAL CROSS-SECTIONS 

The zonally distributed cross-sections of the radiation contributions 
over the ocean-troposphere system are presented in Figs. 18, 19, 20 and 
21 for winter, spring, summer and autumn seasons, respectively. The 
cross-sections, show the results after averaging over the four meridians 
on the mid-seasonal dates conSidered in this study. The results are dis- 
played in the format of Fig. 9. 


In obtaining zonally distributed seasonal means of radiative heating 





rates, denoted Q0(%), at each five-degree multiple of latitude @ in the 
mime 20S,... GSN, values of each radiative parameter listed in Fig. 9 
were arithmetically By areca hover the four meridians. At ¢d = 65N there 
was Only one contribution to the zonal average, while in the Southern Hemi- 
sphere latitudes 20S, 15S and 10S, only two values (on A = 125W and A = 
35W) of each parameter contributed to the means. At 5S, there were three 
meridional sets of radiative parameters entering into Q(¢). Otherwise, 
there were four seasonal values of each radiative parameter entering into 
the arithmetic mean O(¢) at each latitude $? of the cross-sections, Figs. 
ee ep ols 

Consequently, near the northern and southern boundaries of Figs. 18, 
oee,21, the listed seasonal values may not be as representative as in 
mid-latitudes. Nevertheless, the general equatorial-to-polar trend in 


the radiative-change terms appears reliable in both hemispheres. 
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Figure 18(a). 
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Zonally-averaged radiational cross-section for tropical latitudes. 
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Figure 18(b). 


key. 


Zonally-averaged radiational cross-section for higher latitudes. 


All values listed are daily averages in ly min7+ 


Refer to Fig. 9 for 


and are from data for 16 January 1974. 
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Zonally-averaged radiational cross-section for higher latitudes. 
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Refer to Fig. 9 for 


-l and are from data for 16 October 1973. 


-section for higher latitudes. 


Zonally-averaged radiational cross 
All values listed are daily averages in ly min 


Figure 21 (b), 


key. 





A zonal cross-section depicting the mean annual radiative distribu- 
tion has been constructed, Fig. 22, by arithmetic averaging over the 
four mid-seasonal sets of results (Figs. 18,...,21). 

It has also been convenient in earlier sections (cf. Tables III, IV, 
VIII) to compute "weighted averages" with respect to latitudes in the 
Northern Hemisphere. Thus a definition of a weighted-average parameter 
which takes into account the number of observations kK. available at each 


latitude has been given (Meyers, 1975) as 





14 4 
ye an (2 253) erk 

_ 1=1 4 j=1 kK. 1 

So = 16 ee a (6-1) 

Weeeavd. 14 
ye 
iz. 4 OOS 95 
Here etl is the Q-value on meridian j at latitude . and Ks = 1, sce es 


the number of meridional obServations available for the arithmetic 
average Q(0.)-. Note that i = 1,...,14 corresponds to the 14 latitudes, 


>. = Oe ON « 


B. ANNUAL RADIATIVE BALANCE FOR THE EARTH-TROPOSPHERE SYSTEM 

For purposes of this summary, the term BALT (at k=2) in Fig. 9 has 
been redesignated for simplicity as Ri (9), and annual values of ES have 
been plotted in Fig. 23 as a function of latitude. R, represents the 
Deemed Meee clualister Late across the tropopause, that is at the top of 
the troposphere-ocean system considered here. Similarly, the net flux 


at the surface previously denoted BALB in the key, Fig. 9, is 
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Figure 22(a). 


Refer to Fig. 9 


Zonal-annual radiational cross-section for tropical latitudes. 
All values listed are daily averages in ly min7~1 and are computed by arithmetic averaging 


over the four midseasonal results (Figs. 18,...,21l). 
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Figure 22(b). 


key. 


Reter Le Fig. 9 for 
ithmetic averaging 


Zonal-annual radiational cross-section for higher latitudes. 


All values listed are daily averages in ly min7! and are computed by ar 


Over the four midseasonal results (Figs. 18,...,21). 





a SSe_“ 


redesignated here following Malkus (1962) as R($) and its annual distri- 
bution has been graphed against latitude in Fig. 23. Finally, the rela- 


tionship between R. and R in any column is given by 


R =R+R (6-2) 


where R, is the overall net cooling rate or flux-divergence of the tropo- 


eemeric column, e.g., 
R. = BAL24 + BAL46 + BAL6S + BAL810 (6-3) 


foerig., 23, values of R_ (9) have also been plotted against latitude; how- 


ever, a simpler method of obtaining Ro from (6-2), namely, 
R =R -R (6-4) 


has been employed. The resulting zonal annual distributions of each of 
the three parameters, Rue R and Rv: are shown as functions of latitude 
in Fig. 23, where they are superimposed against similar functions drawn 
Peom Sellers (1965). 

The comparison which follows focuses on the annual distribution ie 
the Northern Hemisphere only. It should also be noted that Sellers' 
radiative parameters were taken from climatology at the top of the mean 
zonal atmosphere, whereas this radiative model gives corresponding results 
for level k = 2, based upon tropospheric soundings over an oceanic surface 
Smiley. 

Throughout the latitude range (0-65N), the radiative net flux R at 


the ocean surface is substantially greater than the climatological amount 
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Figure 23. Radiative net fluxes at the tropopause (R,) at the ocean 
Sieeace kjeand net-flux divergence for the tropospheric column {(R_) 
Peonmenewzenalvannial results (Hig. 22(a,b))- Solid lines indicate 
values of Ree R and R, aoeemwesellers “loose 
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a 


of Sellers. This is attributable to the reduced model-albedo values of 
Sec. IV.F. This in turn allowed more insolational-absorption in the 
ocean thereby increasing R over all latitudes as compared to Sellers' 
values (Fig. 23). The term R-R, (which equals -R,) may be regarded as 

an atmospheric "greenhouse effect" in contributing to warming of the sur- 
face both for the model and for Sellers' climatology. For each latitude, 
Fig. 23 indicates that the “model greenhouse" is in close agreement with 
that of Sellers. 

The R_-distributions of both Sellers and of the present model are in 
good agreement across the entire range of Fig. 23, both representing cool- 
ing rates in the troposphere for all ¢. A minimum value of Ro at latitude 
20N appears on both curves and is presumably due to a local maximum of 


IR net-flux divergence F* = Filo’ (resulting from the minimum opaque 


eloud covem in subtropical latitudes. From 50N to 65N, an increase in 


the model] values of Ro appears in Fig. 23 and seems to indicate a reduced 
IR net-flux loss associated with high CL-values (Table III) at these 
northerly latitudes. Sellers" results indicate the opposite trend, i1.e., 
increasingly negative R_-values from 50-60N, but his climatology may not 
include the kind of detail that would indicate reasonable oceanic cloudi- 
ness over these latitudes. 

The radiative model gives R,-values which are generally larger, both 
over tropical latitudes as well as over the northern latitudes, than the 
R, -values of Sellers. As in the R(d) comparison this is a result primar- 
ily of the reduced model-albedo. In the mid-latitude zone 25-40N of 
Fig. 23, close agreement occurs between the two R -curves presumably due 


to such similar climatological effects as for example, the mean polar 
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front. The poSitive-to-negative Fe crossover-point occurs very near to 
40N for both R,-distributions of Fig. 23. The model weighted-average 1s 
R, (6) = 0.048 ly ee if averaging iS considered only to 65N. If nega- 
tive Ry annualovalues Cxist in latitudes © > 65N, an ocean-troposphere 
radiative balance should very nearly be established over the Northern 
Hemisphere. However data for the model were not analyzed over ice-covered 
regions poleward of 65N, where presumably higher surface albedos and 
cooler drier soundings would enable Re to assume larger negative values. 
The radiative-model values of R, (>), R(d) and R_ (9) are in agreement 
with present albedo and IR net-flux observations from recent satellite 
climatology. It should be noted that Sellers' radiative results are 
based on climatology of London (1957), and include also data from the 
land areas of the Northern Hemisphere; on the other hand the radiative 


model values have been based only upon soundings over the ocean. 


SC ROSO- SEASONAL EFFECTS IN THE NORTHERN HEMISPHERE 
ee ON mae Nope lux Bcross the TrOpe pause, Ry 

The seasonal distribution of the tropospheric net flux is pre- 
sented as a function of latitude in Table XI(a) where the listed R, ~values 
correspond to the tropopause, level k = 2. One may discern two geographic 
regimes of annual waves in Table XI(a); a polar zone (25-65N) anda 
tropical-subtropical zone (0-25N). 

The cross-Seasonal variation of Re in the northern regime has 
the following property: R, values are a Minimum in winter increasing to 
a maximum in summer followed by a decrease in autumn. In the tropical 
zone (O-10N), R, -values have positive peaks in spring and autumn. The 


summer minimum of Ry in O-10N is presumably attributable to increased 
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eloudiness associated with the ITCZ in this latitude zone. The subtropi-~ 
cal region 15-25N does not have a minimum in midsummer, and may be simply 
classified as intermediate in behavior between that exhibited by R, in 
the polar and the tropical regimes. These geographic cross-sectional 
effects (Table XI(a)) tend to confirm the general validity of the radia- 
tive model. 
Bec elie Sea-surtace Model) Balance, RK 

The seasonal distribution of the net radiative flux R at the 
Seer s Surface 1S listed as a function of latitude in Table XI(b). Again 
a@ polar zone (25-65N) and a tropical-subtropical zone (0-25N) are defined 
for discussing cross-seasonal effects on R. 

The model results for R are analogous to those just specified for 


R namely: 


ral 
(i) in the polar latitudes there is a single sine wave with 
annual periodicity and maximum in midsummer; 
Gly) ete ron calm lattcudes cnn iene Gcouble peaked character 
with maximum in spring and autumn appears; 


(iii) there exists an intermediate subtropical zone of R-values 


which no longer exhibits a minimum in summer. 


D. TOP OF THE ATMOSPHERE COMPARISON OF MODE] WITH SATELLITE DATA 


i ceerdcadtive Transter Race aeelep of the Atmosphere 


The three parameters that describe the radiative flux at level 


k = O are associated in the following way: 
ee = Oem a : (6-5) 


Qn the daily averaged incoming solar flux, can be related to QAVE and 
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OREF (Fig. 9) by the relationship 


2200 00 
Oe 2) Ale) GAVE ~ OREP (626) 





where QAVE and QREF have already been evaluated according to Eqs. (5-3) 
and (5-6) respectively and the factor 2.00/1.92 has the effect of deriv- 
ing the mean daily solar insolation at level k = O before ozone-oxygen 
absorption. Po the IR net flux to space, was previously discussed in 
cece tits. [Ea. (3-11)])>] and Rays the difference between Qn and Po 
represents the daily-averaged radiative net flux at the top of the 
atmosphere. In order to more easily distinguish between model-values 
Gr Ry and the same parameter from satellite climatology of Raschke et al. 
(1973), Ry is further identified as R MOD er R RAS, respectively. OQOAVE 
is assumed identical both for the mcdel and for the treatment of the 
Raschke data in this study. 

2. Zonally-Averaged Computations of Ry 

With the use of Eqs. (6-6) and (6-5), values of O. ae and Ry 
for each of the four mid-seasonal dates and by five-degree latitude 
intervals have been computed both for the model-parameters and the satel- 
lite-observed (RAS) parameters. The results of the computations both 
seasonally and for the annual mean casSe are presented in Table XII. 

A comparison of Ry MOD with R RAS from the zonal-annual results 
at the bottom of Table XII shows that the crossover from positive-to- 
negative values of Ry occurs near 40N for both model- and Raschke-results. 

Pestaei1stical accuracy check of R MOD as against R RAS data is 
conducted using the standard deviations of the difference parameter 


(R. MOD - R RAS) over all Northern Hemisphere latitudes. A summary of 
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At each latitude the 


“l) for both the model- and Raschke-calculat 


Zonally averaged net-flux parameters at the top of the atmosphere. 


TABLE XII. 
radiat 


LOnS 


(e) denotes annual 


zi @| 


tes the d 


are listed (ly m 


tribution over seasons while 


LS 


Mnadica 


(da) 


‘ 


v 


6) 


"7 


td 


dy 


5 


ive parameters 0Q 


( 


The letter code (a), 


Paticude. 


£ 


distribution of the same parameters as a function o 








the mean and standard deviations of this parameter over the latitude 


range of these five cross-sections (Table XII) follows: 


16 Jan. 16 Apr. Jers ie TovOckr. Annual 
RMOD-R, RAS .004 .024 .019 sole ul) ils wee 
Std. Dev. .028 w087 .053 .039 .039 ly ae 


The zonally-averaged mean (R,MOD-R, RAS) 1S positive for each 
season as well as for the annual distribution (whose mean is 0.015 ly 
pein) - This small positive difference was anticipated from considera- 
tion of Sec. IV.F. (Table IX) where the middle-to-high latitude albedos 
were tuned slightly too small. 

The overall standard deviation in (RMOD-R RAS ) or 0-039 Jy, wee 
is of the same order of magnitude as that of QO, MOD-O, RAS . This compara- 
tively small difference in the case of model versus satellite climatology 


Was a considerable improvement over (A,B,C,D) where untuned albedos were 


used. 


Be) 2 LRATOSPHERIC MODEL RADIATIVE BALANCE 
It is clear that the flux-convergence in the stratosphere is given 


by 


eee 


where Ry is given by Eq. (6-5). Note that Re is identical to BALT of 


Fig. 9 which reduces to 


R, = OAVE - QREF - F.* . (6-7) 


Pel 





Forming the difference, a Ree leads to 


ees ee Ve (PS F*) - (6-8) 
The stratospheric flux-convergence, RRs is displayed in zonal- 
annual format in Table XIII. Based on the Northern Hemisphere weighted 
results (bottom line), it appears that a positive net-flux convergence 
exists between O-65N. However the trend in ROR, at higher latitudes 
1s towards increasingly negative radiative values in the stratosphere 
poleward of 65N such that an annual mean radiative balance may be in- 
ferred. For instance, assuming that in the latitude range 65-90N, an 
average zonal-annual value of ROR, Se OI02 Ly ain SxEStSeats7 7 25N, 
then, the cosine-weighted mean is -0.0126 ly ese in the region pole- 
ward of 65N. Thus the summed Northern Hemisphere weighted mean would be 


zero and a radiative balance would exist in the stratosphere. 


F. ZONAL-ANNUAL NORTHERN HEMISPHERE HEAT BUDGET 
1. Tropospheric Radiation Budget 
The radiative heating rate of a tropospheric column may be ex- 


pressed as a function of R_ (9) by the right side of Eq. (6-9) 


0 OS code 5 ee) 
a a 


) (ea) 
Va 


. 


a is the storage heating rate of the column and oie is the required 
flux-convergence of Sensible and latent heat energies compatible with 
the heat balance at latitude @. R. is the tropospheric radiative net 


Sootmgqucatce ama 1S depicted ian Figs 23. (EB + H.) includes the latent 


ik 
and sensible heat parameters respectively, but these parameters are 


not part of this study and will be considered to be contained in oe 


Simeecimo—-) stor Simplicity. 


Peed 





leat. R R Ro = Re 


N te 

20S eo « LO30 20149 
i> TZ60 SOS) Os 
10 Si 293 -1074 MOAI, 

5 ~ 1456 ©1248 20208 

0 ~1443 wich 70226 

5 Beles) LOY ~0240 
10 e250 ONS S, OZ or 
1s) - 1084 20326 O25 
20 O07 0644 O25 0 
2D -0603 s0367 202316 
a0 ~0441 20293 -0148 
35 S022 Omens Poles 
40 -0031 20052 0029 
45 = (OS =-.0107 -.0030 
50 SOS: = OLE = — Ole 
Do =.0426 = Oseueys: ay Ol 8) 
60 = Volos: = OASIe = VAS 
65N -.0563 = ULOG = ag 
ING. AVG. 70602 -0476 2OlL2G 


TABLE XIII. Computed radiative net fluxes at the top of the 
atmosphere (R_) at the tropopause (R,); and the 
stratospheric absorption (R_-R_). atl parameters 
are listed in the zonal-annual format in ly min7l. 
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The Northern Hemisphere cosine-weighted R, turns out to be 
— y —] 
R, = - 0.1551 ly min 


and if it is temporarily attributable to mean storage-cooling S_ of the 
troposphere alone, the mean storage rate corresponds to a temperature- 
Change rate given by 


1440(R ) 
AP ab 


or ersten || 


ve °C day (6-10) 


with AP ib = 800 mb in the troposphere. The resultant cooling rate (6-10) 
Over the tropospheric depth, 800 mb, with zero lateral flux divergence, 
1s approximately -1.14°C per day averaged over the mean ee tropospheric 
eolumn. 

Thus the general circulation of the atmosphere would have to 
function to bring about the atmospheric heat balance. This process 
would require the proper flux-convergence oa to offset the annual radia- 
tive loss of the average saa column of the troposphere. 

2. Zonal-annual Heat Budget of the Ocean 


An equation similar to (6-9) may be written for the zonal-annual 


heat budget for the ocean as a function of R(d). 


Dn By eR Er) (6211) 


Here, So is the storage rate of the ocean water mass and Oke: is the 
required mean oceanic heat flux divergence for a balance at latitude 9. 
The primary ocean mass heating parameter is the net radiative heating 
minx ndt Ene Ocean SUurtace, R (Fig. 23), and(E + Hp) is for convenience 


considered part of Qo0° 


LS 





The Northern Hemisphere annual weighted mean value of R is 
~ _ al 
R = 0.2027 ly min 


which corresponds (Table XI) to a mean radiative heating rate in the 
water-mass column. It should be noted that this mean heating rate is 
noticeably larger than R, and of opposite sign to the corresponding 
tropospheric radiative cooling effect (R. = OS Suey min ~). How- 
ever, if the three-dimensional oceanic flux divergence of both 


sensible and latent heat is considered, the oceanic heat balance should 


result. 
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a a 


Vite CONCLUSIONS 


The major change in this radiative model from the model evolving 
from (A,B,C,D) was the introduction of tuned cloud reflectances. This 
innovation brought the resultant model-albedos into agreement with 
Satellite climatology for comparable data dates. 

Gridpoint comparisons at the top of the atmosphere of model net flux 
Be and of the net incoming model insolation or were made with similar 
parameters from Raschkes!' satellite climatology (1973). Closer agree- 
ment was obtained between Q, MOD and ORAS than in any of the preceding 
Studies (A,B,C,D). Similarly, close agreement in a least-squares 
sense was achieved here between F MOD and PARAS primarily because of 
an improved formulation for the F_,computation than was used in 


earlier studies. Finally the ocean-atmosphere system net flux 


has only a small least-square error when model and satellite results 
are compared with properly stratified gridpoint and time cross-sections. 
In this study, mid-latitude cloud reflectances were also tuned to 
obtain closer agreement between the global albedo and corresponding 
values from satellite climatology. The result was to systematically in- 
crease O., relative to Raschke's values. This result was undoubtedly 
due to a mismatch between the model clouds and those involved in the 
climatological average. The conclusion reached is that if tuning of 


cloud reflectances of this nature is to be attempted in mid-latitudes 


eS 





using FNWC soundings, more accurate cloud-parameterization formulas 
should be established against observed cloud-cover amounts for identical 
experiment-periods similar to GATE. If then further cloud-reflectance 
tuning 1S necessary, it could proceed from a more certain knowledge of 


the “ground truth” provided by these experiments. 
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COMPUTER PROGRAM 


APPENDIX A. 
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